Valley symmetry breaking in bilayer graphene: a test to the minimal model 
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Physical properties reflecting valley asymmetry of Landau levels in a biased bilayer graphene under 
magnetic field are discussed. Within the 4— band continuum model with Hartree-corrected self- 
consistent gap and finite damping factor we predict the appearance of anomalous steps in quantized 
Hall conductivity due to the degeneracy lifting of Landau levels. Moreover, the valley symmetry 
breaking effect appears as a non-semiclassical de Haas- van Alphen effect where the reduction of the 
oscillation period to half cannot be accounted for through quasi-classical quantization of the orbits 
in reciprocal space, still valley degenerate. 

PACS numbers: 73.43.Cd, 71.70.Di, 81.05.Uw, 72.80.Le 



Graphene, a two-dimensional hexagonal crystal of car- 
bon atoms, has attracted enormous attention in recent 
years [![. Its quasiparticles are massless Dirac fcrmions 
propagating with a velocity 1/300 of the speed of light. 
This allows us to discuss and even observe a variety of pe- 
culiar (ultrarelativistic) phenomena in a condensed mat- 
ter system, such as the anomalous quantum Hall effect 
(QHE), the universal minimum conductivity, Klein tun- 
neling, Zitterbewegung and Schwinger's pair production 
1- 

Departures from strict two dimensionality promises in- 
teresting physics as well. A bilayer graphene (BLG), 
which consists of a pair of single graphene sheets bound 
by weak interlayer Van der Waals forces, has also been 
studied extensively 0, 0] • Its low energy excitations pos- 
sess chiral symmetry and a quadratic spectrum, thus 
combining Dirac and Schrddinger like features. More- 
over, BLG offers interesting opportunities regarding de- 
vice applications. Breaking the layer symmetry opens 
an energy gap between valence and conduction bands 
which in addition can be tuned by electric field effect 
|H, 0, 0, 0] — a clear advantage compared to current 
technology semiconductors. 

Establishing on a firm ground what is the minimal 
model describing the system is an important step towards 
real device applications. For instance, in the early exper- 
iments in BLG the gap was shown to be tunable between 
zero (gapless) and mid- infrared energies [B|, @] , which has 
been confirmed very recently by infrared spectroscopy 
Ell- The relevant minimal model has thus to cor- 
rectly describe BLG up to mid-infrared energies, which 
is a stringent restriction to the available possibilities. 

In the present paper we address the behavior of BLG 
in high magnetic and electric fields applied perpendicu- 
larly. Using as a minimal description the 4— band contin- 
uum model with Hartree-corrected self-consistent gap we 
predict new physics to take place at high enough fields. 
Specifically, for magnetic fields between ~ 20 — 30 T and 



electric fields between ~ 0.5 — 1 Vnm 1 (i.e. densities 
- 5 - 10 x 10 12 cur 2 in the standard 300 nm Si0 2 back 
gate setup), both within experimental reach, we have 
found that a significant valley asymmetry is present with 
the following consequences: (i) the period of de Haas - 
van Alphen (dHvA) oscillations halves, in a clear non- 
semiclassical behavior; (ii) the QHE shows a new quanti- 

2 

zation rule for the Hall conductivity given by a xy = 2^-n, 
with n — 0, ±1, ±2, ±3 . . . Observing experimentally this 
new behavior would put on a firm ground our current un- 
derstanding of BLG. 

Model. — We adopt the 4— band continuum model as a 
key ingredient in our minimal description. The Hamilto- 
nian for one of the two valleys (K point) is given by the 
following 4x4 matrix 0, 0, dH whose elements corre- 
spond to the A and B sublatticcs of the top and bottom 
layers (A lt B x , A 2 , B 2 ), 
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where tt± = Tr x ±m y and the vector tt = p+eA is the mo- 
mentum operator in a magnetic field V X A = (0,0, B). 
The parameter v ~ 10 6 ms _1 is the single layer Fermi 
velocity and t± ~ 0.3 eV the interlayer hopping energy 
0- Note that the zero entries in Eq. fl} can be filled 
with next-nearest-interlayer hoppings originating trigo- 
nal warping and electron-hole symmetry breaking effects 
1 ll | . At the energy scales we are interested in here, how- 
ever, these terms can be neglected in a minimal descrip- 
tion 0, Hill- The parameter A accounts for the layer 
asymmetry induced by an external perpendicular elec- 
tric field E = (0,0, £ e xt)- One way to relate A and the 
applied electric field £ oxt is by writing the local poten- 
tial as a sum of two opposite contributions, one coming 
from the electrostatic energy due to £ e xt which tries to 
polarize the system and a counteractive one originating 
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in the screening properties of the system. In the Hartree 
approximation A can then be written as @, 0? ] 



2 A = e£ cxt c + 6 C ° Art, 

ZEoSr 
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where cq — 3.4 A is the interlayer distance, Eq is the per- 
mittivityof vacuum, e r is the relative permittivity of the 
system [13J , and An = n top — ribottom is the charge carrier 
imbalance between top and the bottom layers. A self- 
consistent procedure is then followed since An depends 
directly on the weight of the wave functions in each layer. 
Instead of £ e xt we use the density n as an externally tun- 
able parameter. The two are related as £p Yt , = en/ (2e £v) 
in the standard back-gate configuration [14|. 

In zero magnetic field the system has a dispersion re- 
lation with the characteristic double-minimum structure 
[l8| . In the presence of a finite magnetic field we can di- 
agonalize the problem by going to the Landau level (LL) 
basis Using the Landau gauge A = (—yB, 0,0), 

and noting that the commutation relation between the 
momentum operators is [^±,71"^] = ^f2eBh, we define 
creation and annihilation operators of the harmonic os- 
cillator as 7T± — ► \/2f ^ and tt t \/2f a for eB ^ 0, 
where I 2 = H/\eB\ is the magnetic length. Eigenvalues 
and eigenstates of Eq. (fTJ) are then obtained by assum- 
ing the wave function to be a linear combination of the 
number states of the harmonic oscillator |m), with inte- 
ger m > 0. The LLs , labeled by an integer k > 0, are 
given by 



V2hv 



H = (si,s 2 ), 



(3) 



with r = -J — 1± and d = —i — A. The label a = (s\, s?) 
specifies the outer and the inner bands (s\ = ±1) and 
positive and negative (s 2 = ±1) energies, respectively. 
For k > 2 the parameter A is given by the roots of the 
fourth-order polynomial 

A 4 - (2k - 1 + r 2 + 2d 2 )A 2 - 2dA 

+ [fc(fc-l) + d 2 (d 2 +r 2 -2fc + l)] =0, (4) 

and the Landau states are 



\k,ri)=[a£\k-l), ft\k), l£\k-2), 6$\k-l) 
where the coefficients are given by 0% = x ~!f^ d a^, 7^ 
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{\l+d)2-{k-i) a k-> °k - (\£+dy--(k-i) a k ana me normal 
ization condition. In addition there are the LLs k = 1 
and k = 0. In the former case A is given by the roots of 



A 3 -dA 2 + (d 2 +r 2 )(d-A) = 0, 
and the Landau state can be expressed as 



(6) 
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FIG. 1: (color online) . Asymmetric structure of Landau levels 
at K and K 1 points in bilayer graphene with energy difference 
A. 



Note that Eq. ((6]) has only three roots, while /i provides 
four labellings. We reserve the label \i = (— , +) for the 
k = LL, where A = +d provides the only eigenvalue. 
We then write the respective Landau state as |0, — , +)) = 
[0, |0), 0, ] T . The two states and |0,-,+)) 

are degenerate in the gapless case d = Q. For the 
K' point the Hamiltonian is given by Eq. (fT]) with the 
replacement tt± — > 7r T . Note that Eqs. (|4]) and ([6]), which 
determine the eigenvalues, include a linear term in d so 
that LLs for K and K' points are not identical for A ^ 
(see Fig. [1]). The eigenvalues and eigenstates are then 
obtained by replacing d — > — d in Eqs. (T5|) and © and by 
reversing the order of the elements in vectors © and ijTj). 
As for the k — LL at K' we reserve the label /1 = (— , — ) 
since the respective eigenvalue has A = — d, and write the 
eigenstate as |0, -, -)) K . = [ 0, 0, |0), 0] T . 

The final ingredient in the present model is the LL 
broadening. This is taken into account through a con- 
stant imaginary retarded self-energy E rot = —iT. Re- 
garding the relevant value of T we note that the dominant 
source of scattering originates from ripples and charged 
impurities (20| . Since BLG is less rippled than its mono- 
layer counterpart, and charged impurities do not pene- 
trate in between the layers, they mainly affect the layer 
closer to the substrate, we take T ~ 0.01ij_ as a conserva- 
tive estimate. This phenomenological value is to be ex- 
pected for a system ballistic on the submicrometer scale 
and, moreover, when used to get the longitudinal resis- 
tivity (by inverting the conductivity tensorgiven below) 
it compares well with experiments in Ref. [3j. 

Self- consistent gap. — The gap A defined in Eq. (|2|) is 
obtained by writing the charge imbalance An in terms of 
the LL weight in different layers, 
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where the first and the second terms denote contributions 
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FIG. 2: (color online), (a) A vs n for different values of 
magnetic field. The inset shows the case of a fixed top-gate, 
(b) Magnetization vs n for different values of magnetic field. 
The vertical lines signal the semi-classical period. 



from K and K' points, respectively. The factor \/{irl 2 ) is 
the degeneracy of LLs per system volume. The function 
fr(E) is the Fermi distribution function in the presence 
of LL broadening, 



fr(x + n) ~ -- -Tan *~ 

Z 7T 1 
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valid in the low temperature limit k B T <C T. In Eq. (|9|) 
the chemical potential /i is given by 



oo 
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(10) 



for a given density n, where fr(x) = fr(x) — fr(x + /i). 
The solution is obtained numerically b y so lving Eqs. ([2]) 
to pop till self-consistency is achieved [iff. 

In Fig. [2ja) we show the self-consistent asymmetry gap 
A at T = as a function of electron density n. Increas- 
ing the applied magnetic field increases screening around 
n w 0, which can be traced back to the valley (layer) 
asymmetry of the th and 1 st LLs [see Eq. ([7])], In par- 
ticular, the th LL appears fully polarized (layer 1 at K 
and layer 2 at K') thus giving the dominant contribution 
to screening whenever the Fermi level lies in between the 
two. The inset shows the case of a finite top-gate voltage 
fixed at n tg — 2 x 10 12 cm 2 The case of fixed fill- 

ing factor and varying magnetic field was considered in 
Ref. [ItJ without LL broadening. 

dHvA effect. — As a consequence of LL quantization 
the magnetization due to the orbital motion of electrons 
shows periodic oscillations as a function of inverse mag- 
netic field - dHvA effect. Equivalently, in systems where 
the carrier density can be tuned, similar oscillations show 
up as /i changes. Here we address the dHvA effect in BLG 
at fixed magnetic field B and varying n through a back- 



gate voltage, taking into account that £ e xt changes as the 
back-gate voltage is tuned. 

The magnetization for fixed density n is given as AI = 
—dil(B, fi(n))/dB, where is the Gibbs free energy and 
fi(n) is the chemical potential satisfying Eq. (fTT))) [ill ]. 
The Gibbs free energy in the presence of LL broadening 
is obtained after standard treatment and is given by 
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(11) 

where f(x) = (c^ x + 1) 1 is the Fermi distribution func- 
tion with /3 _1 = k B T and g(x) = it — 2tan _1 (x/r), the 
volume (area) of the system is V, and v includes not only 
the band indices (si, s 2 ) but also the valley index. With 
the experimental conditions T 3> k B T in mind we con- 
sider the T — version of Eq. (fTl"]) . which still has to 
be evaluated numerically along with the magnetization. 
In particular, summation over the LLs is performed by 
introducing a finite LL cut off k c . We choose discrete val- 
ues of the magnetic field B so that the number of states 
taken into account to obtain 51 is preserved. Since k > 1 
LLs have the same degeneracy proportional to the mag- 
netic field, and since the degeneracy of k = and k = 1 
LLs add up to the same value, it is enough to choose the 
magnetic field B^ a discrete values as k c Bk c = Const [53 ] . 
A similar idea has been used for monolayer graphene in 
Ref. 0. 

In Fig.^b) we show the dHvA effect as the carrier den- 
sity is changed at magnetic fields B = 10, 20,30 T. The 
vertical dashed lines signal the oscillation period given by 
the semiclassical approach. The latter applied to a two- 
dimensional system in constant magnetic field and vary- 
ing density implies: a peak in magnetization whenever 
the semiclassical orbit coincides with the Fermi surface; 
a complete oscillation whenever the area enclosed by the 
Fermi surface A changes by the quantum of area enclosed 
by the semiclassical orbit, AA = eB/h pil. Eij]. Owing 
tho the cylindrical symmetry of the bands [Eq. ([I])] the 
oscillation period can be expressed in terms of carrier 
density as 



(12) 



where we used A = |?7|/(47r) and included spin and 
valley degeneracies. When n (i.e. the gap) is small 
the oscillations in magnetization follow the semiclassi- 
cal prediction [Fig. [2jb)] . With increasing n (i.e. larger 
gap) we observe deviations from the semiclassical value 
[arrows in Fig. H{b)] and the new oscillation period is 
half of Eq. (fT2j) . This provides a clear sign of the val- 
ley symmetry breaking, whose experimental observation 
could be possible by diamagnetic current measurements. 
Moreover, oscillations with a period which is half of 
Eq. (fT2"j) should also be observed in longitudinal conduc- 
tivity (Shubnikov-de Haas oscillations). 
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FIG. 3: (color online). Hall conductivity vs electron density 
at different magnetic fields: B = 10 T (a), 20 T (b), 30 T (c). 



QHE. — The QHE is discussed by extending the results 
of Ref. [26| for the gapless BLG case. We use the Kubo 
formula to get the longitudinal and Hall conductivities, 



self-consistent gap and Landau-level broadening we have 
shown that the effects of valley symmetry breaking are 
manifested in non-scmiclassical oscillations of the mag- 
netization and anomalous quantized Hall steps. Their 
experimental detections is within reach, and would allow 
for a critical test of the minimal model for BLG. 

We thank M. Koshino and M. A. H. Vozmediano for il- 
luminating discussions. MN is supported by Global Cen- 
ter of Excellence Program "Nanoscience and Quantum 
Physics" . EVC is financially supported by the Juan de la 
Cierva Program (MCI, Spain). BD was supported by the 
Hungarian Scientific Research Fund under grant num- 
ber K72613 and by the Bolyai program of the Hungarian 
Academy of Sciences. MN and EVC acknowledge the 
visitors program at the Max-Planck-Institut fur Physik 
komplcxcr Systcmc, Dresden, Germany. 



Re dij (uj) = 



Im n, j (u> + iry) 



(13) 



where ny(w) = Ily(w) — Hij(0) with {i,j} e {x,y}. 
The polarization function Ey(w) is given by the current- 
current correlation function, and obtained as the analyt- 
ical continuation of the Matsubara form: 
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where /i and v include not only the band indices but also 
the valley index. The matrix 7^ is defined by 7 = V p 7i, 
and HEl = E% with LL Eg as in Eq. ©. We have defined 
\u) n = iu n + [/x + isgn(u; n )r]//i and w+ = w n + v m with ui n 
(v m ) as the Matsubara frequency of fermions (bosons). 

Numerical results for the Hall conductivity as a func- 
tion of carrier density (no top gate is assumed) are shown 
in Fig. [3] for different magnetic fields ^ . For B = 10 T 
the Hall conductivity follows the gapless BLG quantiza- 
tion rule a xy = Ane 2 /h with n = ±1, ±2, ±3 . . . 0). For 
B = 20 — SOT new quantized Hall steps [dashed hori- 
zontal lines in Figs.^b) and[J2c)] appear at high carrier 
densities (large asymmetry gap). The new quantization 
rule reads a xy = 2ne 2 /h and is a direct consequence of 
the degeneracy lifting between LLs from K and K' . In- 
cluding a top gate enables the observation of the new 
quantization rule near and at the neutrality point [HI]. 

Conclusion. — A remarkable property of BLG is the 
possibility to open and tune an energy gap by break- 
ing the layer symmetry. At finite magnetic fields this 
layer asymmetry translates into an asymmetry between 
the two valleys, no longer protected by time-reversal 
symmetry Using the 4-bands continuum model with 



[1 
[2 
[3 
[4 

[5 

[7 

[8 
[9 
[10 

[11 
[12 

[13 
[14 



[15. 
[16 

[17; 

[is: 

[19 

[20 

[21 

[22 
[23 

[24 
[25 

[26 

[27' 



A. K. Geim, Science 324, 1530 (2009). 

A. H. Castro Neto et al, Rev. Mod. Phys. 81, 109 (2009). 

K. S. Novoselov et al, Nat. Phys. 2, 177 (2006). 

E. McCann and V. I. Fal'ko, Phys. Rev. Lett. 96, 086805 

(2006) . 

T. Ohta et. al, Science 313, 951 (2006). 

E. McCann, Phys. Rev. B 74, 161403 (2006). 

E. V. Castro et al, Phys. Rev. Lett. 99, 216802 (2007). 

J. B. Oostinga et al, Nature Mat. 7, 151 (2007). 

Y. Zhang et al, Nature 459, 820 (2009). 

K. F. Mak et al, Phys. Rev. Lett. 102, 256405 (2009). 

J. Nilsson et al, Phys. Rev. B 78, 045405 (2008). 

E. V. Castro et al, arXiv:0807.3348. 

The value of e r depends on the substrate. For a Si02 / air 
interface e r ~ 2.5. Here we use e r ~ 1 as a reference. 
In dual-gate cases £ cx t and n are independent: £ cx t = 
e(n bg — n tg ) / (2eoe r ) and n = n hg + n tg , where n bg and 
ntg are back- and top-gate induced densities pj|. Il5l[. 
S. Kim and E. Tutuc, arXiv:0909.2288. 
The LL sum in Eqs. © and (Tj0| is cutt off at k c w 1000. 
A similar approach was followed in Ref. fil\ . 
M. Mucha-Kruczynski, E. McCann, V. I. Fal'ko, Solid 
State Commun. 149, 1111 (2009). 

F. Guinea, A. H. Castro Neto, N. M. 
Rev. B 73, 245426 (2006). 
J. M. Pereira, Jr., F. M. Peeters, and 
Phys. Rev. B 76, 115419 (2007). 

M. I. Katsnelson and A. K. Geim, Phil. Trans. R. Soc. A 
366, 196 (2008). 

S. G. Sharapov, V. P. Gusynin, and H. Beck, Phys. Rev. 
B 69, 075104 (2004). 

We use k c ~ 1000 and checked consistency of results. 
M. Koshino and T. Ando, Phys. Rev. B 75, 235333 

(2007) . 

L. Onsager, Philos. Mag. 43, 1006 (1952). 

I. M. Lifshitz and A. M. Kosevich, Zh. Eksp. Teor. Fiz. 

29, 730 (1955). 

M. Nakamura, L. Hirasawa, and K. Imura, Phys. Rev. B 
78, 033403 (2008). 

Summing over 100 LLs is enough for the conductivity. 



R. Peres, Phys. 
Vasilopoulos, 



